Transport properties of bottomed mesons in a hot mesonic gas 
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In this work we evaluate the B-meson drag and diffusion coefficients in a hot medium constituted 
of light mesons (tt, K, K and rj). We treat the B-meson and _B*-meson interaction with pseudo- 
Goldstone bosons in chiral perturbation theory at next-to-leading order within the constraints from 
heavy quark symmetry, and employ standard unitarization techniques of NLO amplitudes in order 
to account for dynamically generated resonances (leading to a more efficient heavy-flavor diffusion) 
and thus reach higher temperatures. We estimate individual meson contributions from the gas to 
the transport coefficients and perform a comparison with other findings in literature. We report a 
bottom relaxation length of 83 fm at temperature of 150 MeV and for typical momenta of 1 GeV, 
at which our approach is reliable. Compared to a charm relaxation length of 40 fm in the same 
conditions, we conclude that the B mesons provide a cleaner probe of the early stages of a heavy-ion 
collision. 



I. INTRODUCTION 

The features of matter formed in Heavy Ion Collisions 
(HICs) have been a subject of great interest in the last 
decades. In this scenario, heavy-flavored hadrons play an 
essential role since they carry heavy quarks produced in 
the early stage of the collisions. Therefore, heavy mesons 
are interesting probes to understand the evolution of par- 
tonic matter since its creation, unlike pions and kaons 
which can be produced in the thermal medium at later 
stages. 

However, it is worth noticing that the momentum spec- 
tra of charmed and bottomed mesons extracted from 
HICs undergo modifications due to their interactions 
with the hadron medium, constituted of pions and other 
particles. In this sense, the diffusion of heavy mesons in 
an equilibrium hadronic gas must be taken into account, 
and may be properly studied in the framework of kinetic 
theory to compute transport coefficients. 

Different approaches have been used to study various 
aspects of this topic [JQ- In particular, in Ref. Q heavy 
quark effective theory (HQET) and chiral perturbation 
theory (ChPT) have been employed, focusing on the low- 
est possible temperatures. Attempts to reach higher tem- 
peratures in the context of charmed mesons close to the 
crossover to the quark and gluon plasma have been done 
in Refs. @, Q. In addition, in Ref. @ the transport coef- 
ficients of i?-mesons have been obtained with the use of 
scattering lengths as dynamical input. 

In our recent work 0] the transport coefficients of 
charmed mesons in a hot pion gas were computed, ex- 
ploiting ChPT at next-to-leading order (NLO) and em- 
ploying standard unitarization as guiding principle to 
reach higher temperatures and account for the contri- 
bution of resonant channels. Thus, a natural question 
arises about the application of this approach to bottomed 
mesons, which would allow a comparison with existing 
literature and, hence, a better comprehension of this is- 
sue. 

In the present work we extend the framework used in 



[H to evaluate the i?-meson drag and diffusion coeffi- 
cients in a hot mesonic gas, including pions, kaons and 77 
mesons, with special attention to the contribution of the 
heavier states with respect to the pion gas. We perform a 
detailed analysis of the temperature and momentum de- 
pendence of these coefficients as well as their static limit 
(vanishing heavy-meson momentum) and scaling proper- 
ties with the heavy-meson mass. 

The organization of this paper is as follows. In Sec. [Hi 
we introduce the transport coefficients and the chiral La- 
grangian density that describes the interactions between 
B-mesons and light mesons. Afterwards, Sec. IIIII is de- 
voted to obtain the unitarized scattering amplitude and 
fit the relevant free constants to available data. The 
transport coefficients are evaluated and analyzed in de- 
tail in Sec. IIVI for a pure pion gas. A pertinent discus- 
sion of the role of unitarization of scattering amplitudes 
in heavy-quark diffusion is also contained in this section. 
The modifications in the transport coefficients due to in- 
clusion of kaons and 77 mesons in the thermal bath are 
studied in Sec. |Vj A summary and concluding remarks 
are given in Sec. I VII 



FORMALISM 



II. 



A. Transport coefficients 

We evaluate the transport coefficients of stable B- 
mesons propagating through a hot mesonic gas by using 
the scattering amplitudes obtained from chiral pertur- 
bation theory, taking into account unitarity and heavy- 
quark symmetry. We assume that the density of pseu- 
doscalar (B) and vector (B*) bottomed mesons is very 
small, so we will neglect collisions among bottomed 
mesons themselves and concentrate only on their inter- 
action with the light meson gas in thermal equilibrium. 

Following the assumptions made in Ref. the 
momentum-space distribution of bottomed mesons must 
relax via the Fokker-Planck equation. We choose the mo- 
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mcnta of the clastic collision between a meson and 
a light meson <j> as 

fl«(p) + 0(q) -> B«(P " k) + </>(q + k) . (1) 

In this context, the evolution of the momentum distri- 
bution of bottomed mesons due to their interaction with 
the isotropic mesonic gas is fully controlled by the drag 
(F) and diffusion (Fq,Fi) coefficients, written as 



dk w(p, k) — =- 



Hp 2 ) = 
r ( P 2 ) - j / dk w ( P ,k) 

Ti(p 2 ) = i / dk W (p,k) 



(2) 



where w(p, k) denotes the collision rate for a bottomed 
meson with initial (final) momentum p (p — k) , 



w(p, k) 



9<t> 



U(q) [l + / (q + k)] 



2E$ 2E+ 2E*_ k 2E q+ , 



E 



|X B0 (M,x)| 2 - 
(3) 



In Eq. ([3]), /</>(q) is the bath's distribution function; 
M B4> stands for the Lorentz invariant B meson - light 
meson scattering matrix element, is the Goldstone bo- 
son isospin degeneracy (i.e. = 3 for the pion), and x 
denotes the possible spin degrees of freedom. 

In the next subsection we derive the scattering ampli- 
tude A4 for bottomed mesons in the light meson medium, 
necessary to evaluate the three transport coefficients in- 
troduced above. 



B. Effective Lagrangian for B- meson and light 
meson interaction 



Our task now is to construct the chiral Lagrangian 
density that describes the interactions between the J = 
and J = 1 5-mesons and light mesons. In this sense, we 
note that heavy-meson ChPT, described in Refs. 
for the case of charmed mesons, can be applied to the 
B-meson sector as well. 

Let us start by introducing the pseudoscalar Gold- 
stone bosons. They follow the nonlinear realization of 
the SU{Nj)l x SU(Nf)n chiral symmetry, given in the 
exponential parametrization 



U = exp 



F 



(4) 



with F being the Goldstone boson decay constant in the 
chiral limit and <f> the matrix incorporating the pseu- 



doscalar Goldstone bosons, 



IT 



V 



7T 

K~ 



1 



K° 



K+ 

K° 
2 

V6 



(5) 
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Under chiral symmetry, the matrix U is transformed as 

U ^U' = LUR f , (6) 

where L and R are global transformations under SU (3) l 
and SU(3)n, respectively. The Goldstone boson kinetic 
term of the effective Lagrangian is explicitly invariant 
under this symmetry 

C* = ^{d li tfd»U) . (7) 
For convenience, a matrix u is introduced as 

u = VU , (8) 
and is transformed under chiral symmetry as 

u -> u = LuW^ = WuR^ , (9) 

where W is a unitary matrix expressible as a certain com- 
bination of L, R and <f>. The axial and vector fields are 
constructed as 



u,, 



i (w t 9 A , 



,11 



whose transformation laws read 



r' = WT^ + Wd^W 1 * 



(10) 



111) 



3 



Finally, the covariant derivative which acts on the heavy 
meson field reads 

v M = - r M . (12) 

With these ingredients, the leading-order (LO) chi- 
ral Lagrangian £W involving heavy mesons and pseu- 
doscalar Goldstone bosons is given by 

£(i) = v^PV/ - mlPP 1 - VP*" V^P** 
+m 2 B P*^P*^ + ig (P^u^P 1 - Pu»P* ] ) 

(13) 

where P = (P-,P°,P°) and P* = , P*°, P* )^ 

are the SU(3) anti-triplets of spin-zero and spin-one B- 
mesons with the chiral limit mass tub and ms* , re- 
spectively. This Lagrangian is indeed invariant under 
SU (3) l x SU (3)r symmetry. The NLO chiral Lagrangian 
reads 

= -h a PP\ X +)+hiPx+P ] + /i 2 PP t (^i V ) 
+A 3 Pu'*u M P t + /nV^PV^^tt") 

i/isV/K./jv/ + h p*»pp( x +) 

-hiF*» X +P? - h 2 P*»P*J{u v u v ) 
-hiP^u u u v P^ - /> 4 V M P* Q V„P*t <«"«") 
-AsV^K.u^V^ , (14) 

where 

X+ = ir + M X U (15) 

and x = diag(TO^., m\, 2m|- — mj) being the Goldstone 
boson mass matrix. The twelve parameters hi,hi(i = 
0,...,5) are the low-energy constants (LECs), to be de- 
termined. However, we can make use of some constraints 
to reduce the set of free LECs. First, it should be no- 
ticed that in the limit of large number of colors (N c ) of 
QCD [To| . single- flavor trace interactions are dominant. 
So, we fix ho = h-2 = /14 = ho = h% = /14 = hence- 
forth. Besides, by imposing the heavy-quark symmetry, 
it follows that hi ~ hi. 

In the following, we use the lowest order of the per- 
turbativc expansion of the quantities T^, and x+ m 
Eqs. (TT3"|) and (fl4"|) . and construct the scattering matrix 
of the interaction between the bottomed mesons and the 
pseudoscalar Goldstone bosons. 

III. SCATTERING MATRIX FOR THE 
BOTTOMED MESON IN THE MESON GAS 

A. Scattering matrix elements 

With the Lagrangian in Eqs. (TIB"]) and (TT4")) we are 
able to calculate the scattering of pseudoscalar Goldstone 



bosons <f> off pseudoscalar B mesons as well as vector 
B* mesons. In Fig. [T] we show the tree-level diagrams 
constructed from the LO and NLO interactions. These 
include both contact interactions and Born exchanges. 




FIG. 1: Tree-level diagrams relevant to the scattering ampli- 
tudes for B w </> -»■ processes. The solid, double and 
dashed lines represent the B-mesons, _B*-mesons and Gold- 
stone bosons, respectively. 

Note that the Lagrangian density in Eq. (fT4|) has 
been manifestly constructed maintaining chiral symme- 
try, that is then broken only carefully in perturbation the- 
ory upon expanding in fields and derivatives to construct 
the LO and NLO chiral amplitudes. However, since the 
bottomed mesons are heavy fields, the heavy-meson sym- 
metry should be recovered by taking nij — > 00. We dis- 
cuss below the implications of this limit in the scattering 
amplitude. 

The spin-changing amplitudes B*<j) — > Bcf> and B(j) — > 
B*cf>, shown in diagrams (b) and (c) of Fig. [TJ should 
vanish in the limit — > 00, since a collision with a 
Goldstone boson cannot change the heavy-quark spin, 
that decouples in this limit. It can be easily proved that 
indeed these amplitudes are subleading in 1/ mi, [24| . 

Turning to the elastic B<p and B*(f> amplitudes, dis- 
played in diagrams (a) and (d) of Fig. [TJ one finds [3] 
that the Born exchange terms (proportional to g 2 with 
an intermediate bottomed meson propagator) are sub- 
leading in HQET, and therefore are suppressed by m^ 1 
with respect to the contact interaction. 

Hence, the final amplitude for scattering off a bottom 
quark in the light meson gas, at NLO in the chiral expan- 
sion and LO in the heavy quark expansion, irrespective 
of whether the heavy quark is in a B or a B* meson state, 
is given by 

V 2^(s -u) + y^hi + ^h 3 (p 2 ■ pa) (16) 
2C3 

+ -p2 h 5 [(Pi -P2)(P3 -Pi) + (Pi -Pi){p2 -Pa)] 1 

where C, (i — 0,...,3) arc channel-dependent numerical 
coefficients in isospin basis, as collected in Table U (chan- 
nels are denoted as B^cj)(I) 7 with total isospin /). 
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1/3 
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-1 


1 
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1/3 



TABLE I: Coefficients of the scattering amplitudes for the 
B^(f>(I) channels with total isospin / in Eq. (|16l) . 



B. Unitarized scattering amplitude 



It is well known that ChPT, involving a perturbativc 
expansion up to a certain order, is bound to work prop- 
erly only at low energies. Particularly, it cannot de- 
scribe the presence of resonances in specific scattering 
channels, since a resonance shows up as a pole in the 
S'-matrix complex energy plane. Moreover, the use of 
ChPT amplitudes at a finite order leads to cross sections 
which increase monotonically with energy, eventually vi- 
olating Froissart bounds imposed by the unitarity of the 
S'-matrix. This, in turn, limits the aplicability of the the- 
ory in finite-temperature calculations, since the higher 
the temperature means the more energy available for the 
two-body interaction. 

Imposing unitarity of the scattering amplitudes solves 
these problems and furthermore accounts for resonances 
dynamically generated from the LO amplitudes in those 
channels where the interaction is attractive. In the 
present case, and motivated by our former experience 
in the charm sector we can expect to find reso- 
nances in the SWtt(1/2) and the B^K(0), B^K(0) 
S-w&ve channels, cf. Table U (we briefly review the exper- 
imental situation of the i?-meson excitation spectrum in 
Sec. IIII Cj) . As a resonant interaction implies a more effi- 
cient diffusion (and, thus, shorter thermalization times), 
we believe unitarization of the scattering amplitudes is 
mandatory in this approach. 

Following Refs. HG3, 

using on-shell unitarization via 
the Bethe-Salpeter equation, the unitarized scalar ampli- 
tude T can be written as 



dimensional regularization, this integral reads 

i r _2 



G(s) 



__ ia(M)+ l n _ 



[ln(s — (r 



2.s 

2<Zv^) 



In 



+ ln(s + (m| - m\) + 2q^) 



4 

ln(s - (m| - m\) - 2q^/s) 
-ln(s+(m!-m|)-2gVi)-27ri]} , (19) 



where [i is the regularization energy scale, a(/i) is a sub- 
traction constant which absorbs the scale dependence of 
the integral, and q is the modulus of the light meson's 
three-momentum in the CM frame, 



~/=]/i s ~ ( m <t> + m s) 2 ] [s - (m^ - tob) 2 



2^ 



(20) 



It is worth mentioning that heavy-meson spin sym- 
metry guarantees the same scattering cross section for 
both B(j) or B*4> channel, and no further spin averaging 
is needed. Then, we use in the colision rate defined in 
Eq. ([3]) the scattering matrix element given by 



J2\M B 4s,t, x )\ 2 = \T t 



(21) 



where |Tb</>| 2 is the isospin averaged unitarized ampli- 
tude, namely: 



22) 



with T 1 being derived from Eq. (|17|) in the total isospin 
basis. 

For the sake of comparison with other systems of in- 
terest, we also evaluate the corresponding cross sections 
in the CM frame, 



<7B4>(s) = 



16tts 



\Mi 



(23) 



T(s) 



-V(s) 



1 - V(s) G(s) 



(17) 



where V(s) is the S-wave projection of the scattering 
amplitude in Eq. (fTB"]) . and G(s) stands for the two- 
meson loop integral, 



G(s) = i 



d A q 



(27r) 4 (P - q) 2 - m 2 B +ieq 2 - m 2 , + it 



(18) 

where is the mass of the light meson. Employing 



C. Free Constants 

The only remaining task before proceeding with the 
calculation of the drag and diffusion coefficients is to de- 
termine the free constants of the theory from available 
data. 

At the level of precision that we are working, the pion 
decay constant in the chiral limit can be approximated 
by the physical value, F = 92 MeV. The values of the 
meson physical masses that we use are: mg = 5279 MeV, 
m B , = 5325 MeV, m Bs = 5366 MeV, m B , = 5415 MeV, 
m„ = 138 MeV, m K = 496 MeV and m„ = 548 MeV 
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Let us focus on the scattering of pions off B-mcsons. 
We use the renormalization scale fi = 1 GeV, and the 
scheme is such that the subtraction constant a(/i) = 
—3.47. Details in the choice of these numbers are given 
in Appendix [A] For this choice of parameters we show 
in Fig. [2] the squared amplitude, |T| 2 , for Bit scatter- 
ing with isospin / = 1/2, keeping only the LO (s — u) 
term in the elastic amplitude V. This channel, as ex- 
pected, exhibits a resonant behavior with a peak around 
y/s ~ 5530 MeV, which is in good agreement with the 
former determinations in [13Hl5| . 

Notice that we have only considered the isospin chan- 
nel I = 1/2 in the square amplitude displayed in Fig. [2] 
The exotic 7 = 3/2 channel is non- resonant in this case 
and in all situations presented below. 

In addition, we compute \T\ 2 for B*ir scattering with 
/ = 1/2, just by replacing ms by m* B in V(Btt — > Bit). 
The result is plotted in Fig. [3j where now the peak shows 
up at y/s ~ 5580 MeV, also in agreement with Ref. fl6| . 

Some comments about the theoretical and experimen- 
tal knowledge of the i?-meson spectrum are in order: 
From heavy quark considerations (confirmed by the well 



known charmed meson spectrum) one expects a spin 3/2 
doublet and a spin 1/2 doublet as first excitations. The 
former is composed by the states B\ and B^, both of 
which have been observed as the Bi(5721)° and B* 2 (5747) 
resonances. These states are narrow, the B\ state decay- 
ing mostly into B*tt in D-wave. Therefore, this state does 
not correspond with the resonance that is generated in 
our amplitude. On the other hand, other two states Bq 
and Bi are supposed to exist (in analogy with the Dq 
and Di states for charm) that should be much broader 
and have not been discovered yet. The latter decays into 
B*ir in S'-wave, which would naturally correspond to the 
state that we find in Fig. [3] As its mass is still not deter- 
mined experimentally, we can only guess that it is very 
close to the mass of the other B\ state, appearing at 
5720 MeV. We justify this guess based on the fact that, 
in the charm sector, the masses of the two D\ states are 
nearly the same, and thus we expect the same behavior 
for the B system provided that heavy quark symmetry 
works. Many quarks models also predict a similar mass 
for the two Bi . See Tabic [IT] for a small review. 



Heavy Spin, J" 


D (quark model) 


D (experimental) 


(M,r) MeV 


B (quark model) 


B (experimental) 


(M,r) MeV 


1/2, 0+ 


Do 


L>S(2400) 


2318, 267 


Bo 


? 


? 


1/2, 1+ 


D 1 


Di(2430) 


2427, 384 


Bi 


? 


? 


3/2, 1+ 


D 1 


£>i(2420)° 


2421, 27 


Bi 


Si (5721) 


5723, ? 


3/2, 2+ 


D 2 


L>3(2460) 


2466, 49 


B 2 


B* 2 (5747) 


5743, 23 



TABLE II: List of heavy meson states. Left side: D meson. Right side: B mesons. The two states in bold font are the ones 
that we obtain in our unitarization scheme. Experimental data taken from [l2l |. 



The NLO terms, containing the hi constants, may be 
used to improve our results concerning the position of 
the observed resonances in view of the former discussion. 
hi is fixed by the mass difference between the strange 
and non-strange i?-mesons, as obtained from the chiral 
Lagrangian C^ 2 \ We have 

m B s ~ m B = -4hi(m 2 K - ml), (24) 

for the pseudoscalar B mesons, and 

m%. - 171%, = -4hi(m 2 K - ml), (25) 

for the vector B mesons. Replacing into these equa- 
tions the values of the masses introduced above, we get 
hi = —1.020 and —1.064 for the scalar and vector cases, 
respectively. We shall adopt in our computations an av- 
erage value: hi = —1.042. We note that the value of 
hi does not yield any relevant changes in the squared 
amplitudes, since the corresponding term in the NLO 
amplitude is multiplied by a small ml constant. 



The last free LECs to be estimated are /13 and /15. We 
proceed to fit them by demanding that the B*tt squared 
scattering amplitude peaks at the mass of the Bi(5721)° 
resonance [I(J P ) = 1/2(1+)]; (5723.5 ± 2.0) MeV [3. 
Following the discussion in [8j we can estimate a valid 
range for h$ from our results in the D sector. Taking 
into account that /15 scales as /15 ~ w^ 2 , we use for the 
B sector: 

„ f = k? (^)\ <*, 

leading to a value of /15 = —0.04 GeV -2 . The value of 
/13 does not affect very much the pole position of the 
resonance. 

Choosing a value of /13 = 2.5 we obtain the final re- 
sult in Fig. |U The predicted mass for this state from the 
peak of the amplitude is still about 100 MeV below the 
£?i(5721)° mass. We find, as discussed in Appendix [A] 
that it is not possible to fix simultaneously the low en- 
ergy and subtraction constants -in a natural choice- in 
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FIG. 2: Square amplitude for Bn scattering with isospin I = 
1/2 and at LO in ChPT. 
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FIG. 3: Square amplitude for B*n scattering with isospin 
I — 1/2 at LO in ChPT, computed replacing wis by m B in 
Eq.[T7] 



FIG. 4: Cross section for B*n scattering with isospin / = 1/2 
and at NLO in ChPT, computed replacing tub by m* B , and 
considering (hi, h 3 , h 6 ) = (-1.042,2.5,-0.04 GeV~ 2 ). The 
resonance is to be understood as the broad B\. 



close to the results in [17j of —0.17 fm. 




order to get values nearer the expected £?i(5721) mass. 
Furthermore, our findings are in agreement with those 
in Refs. [ID, Q, HBj- Because of the lack of experimen- 
tal evidence of a B\ state decaying into £?*7r, the most 
reasonable choice is to keep the values of the LECs de- 
termined above. 

For completeness, we also compute Bit scattering at 
NLO in ChPT. This channel is also resonant, with mj, = 
5534 MeV. Again, there is no data yet at our disposal to 
compare with. The total cross section for the Bit channel 
(considering the two isospin channels 1 — 1/2 and 3/2) 
is shown in Fig. [5] 

The 1=1/2 cross section at threshold is 18.1 mbarn, 

l li 

or equivalently an S- wave scattering lenght of a B7r = 
0.38 fm or m^aB,, = 0.26 is found, in total agreement 
with the results of [ll| but larger than the result of [l7| , 

1 /2 

a B \. = 0.25 fm. For the repulsive channel / = 3/2 we 
find a cross section of 1.5 mbarn, or ti^t = —0.11 fm, 



FIG. 5: Total Bn cross section (both 1/2 and 3/2 isospin 
channels) at NLO in ChPT, (hi, h s , h 5 ) = (-1.042, 2.5, -0.04 
GeV" 2 ). 

After the estimation of the free relevant constants for 
the B^<j> scattering amplitude, we are now in position to 
compute the transport coefficients, which are the subject 
of the next section. 



IV. DRAG AND DIFFUSION COEFFICIENTS 

We proceed to calculate the F, To and Ti transport 
coefficients defined in Eq. @, for bottomed mesons in 
a pure pion gas. We expect the contribution of pions to 
be the most relevant one because of their large multiplic- 
ity in comparison to other particles. For completeness, 
in this work we also account for the effect of the other 
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FIG. 6: Momentum-space drag coefficient as function of tem- 
perature for several bottom quark momenta in the pion gas. 



members of the light meson SU(3) octet. A detailed dis- 
cussion is done in Section Ivl 

In Fig. |6]we show the dependence on the temperature 
and i?-meson momentum of the drag coefficient F. Wc 
observe an increase of a factor of about 6-8 in the range 
from T = 100 — 180 MeV, which means that the temper- 
ature evolution of the thermal medium can modify the 
nature of its interaction with the i?-mesons. So, the drag 
in a heavy-ion collision is considerably strengthened in 
the hotter stages, with a significant interaction between 
heavy mesons and thermal medium, and accordingly at 
larger momentum transfers. However, as the tempera- 
ture of thermal bath diminishes, the magnitude of inter- 
action decreases and the bottomed mesons move more 
freely. The drag coefficient exhibits a mild momentum 
dependence of about 10% in the range [0.5, 2.5] GeV, 
whereas more pronounced at low p. From this coefficient 
wc estimate the relaxation length of bottom quarks in the 
hadronic medium. Taking the pion gas at a temperature 
of 150 MeV (at which our approach is reliable), we find 



X B {T = 150 MeV,p 



1 GeV) = I 



1 



0.01 



fm = 100 fm 
(27) 

for bottomed mesons travelling with a typical momentum 
of 1 GeV. This value of Xb is considerably bigger than the 
case of charmed mesons, evaluated in the same approach 
in Ref. [|[ (there A^> ~ 40 fm). Thus, it is reasonable 
to expect that the lifetime of the pion gas (typically 5- 
10 fm) is smaller than the relaxation time of both D 
and B mesons, which means that heavy quarks do not 
completely relax before leaving the hadronic medium. 

Another point worthy of mention is that the drag co- 
efficient for charmed mesons computed in our previous 
work [J is about three times larger than the one for bot- 
tomed mesons at small p. This scaling is the correct one 
when looking at the non-relativistic expression of this co- 



F 



:<tP 



w. 1 
T m B 



(28) 



where a is the total cross section and P the pressure of 
the gas. Note that there is a dependence on the heavy- 
flavor mass in the denominator. This makes the drag 
coefficient smaller for heavier mesons. With respect to 
the D meson system, the drag force is suppresed by a 
factor of ms/rn/) ~ 2.8 in very good agreement with 
what we observe in our results. Note that in the low 
energy limit the cross section does not depend on the 
heavy meson mass. This is not the case at higher tem- 
peratures where Eq. (|28p contains further corrections in 
powers of 1/mg. This scaling is also observed for the 
drag coefficient of b and c quarks in the static limit be- 
yond the critical temperature within the phenomenolog- 
ical approach of Ref. [HI, Ojl, in total accordance with 
our findings in the hadronic phase. We observe, however, 
that this scaling is not maintained in the results of [3, H| . 
There, a similar approach based on heavy-meson ChPT 
amplitudes is used to calculate the transport coefficients 
of D- and B-mesons. However, the authors employ NLO 
perturbative amplitudes in the case of Dtt scattering [3j , 
overestimating the effect due to the high-energy depen- 
dence of their cross sections. This is partly solved in Q 
for the B system, where NLO amplitudes at threshold 
energy (scattering lenghts) are used in the evaluation of 
B-meson transport, thus taming the high-energy behav- 
ior of the amplitudes. In contrast, this approximation 
underestimates diffusion from resonant scattering in the 
B^tt{1/2) channel. This discussion reinforces the role of 
unitarization of low-energy scattering amplitudes to ob- 
tain realistic transport coefficients at high temperatures 
in the heavy-flavor sector. A more detailed discussion of 
the use of perturbative versus unitarized amplitudes is 
carried over at the end of this section. 

In Fig. [7] the temperature dependence of the B meson 
diffusion coefficients Tq and Ti is displayed for several 
momenta. The two coefficients become degenerate in the 
static limit, p — > 0. At p = 0.1 GeV this limit is already 
well reached. Interestingly, one finds that numerically 
this coefficients are very similar to those obtained for the 
D meson case in Ref. [4[ • Looking at the non-relativistic 
expression of these coefficients gives a clue for this fact: 



To,Li ~ ^aP\pm^T 



where we have used the Einstein relation 



F 



(29) 



(30) 



to obtain this equation. The coefficients r and Ti are 
independent on the heavy quark mass in the nonrelativis- 
tic limit. For this reason the results in Fig. [7] are basically 
the same as those in Ref. In Ref. OH3 very sim- 
ilar diffusion coefficients arc found in the c and b quark 
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FIG. 7: Momentum-space diffusion coefficients as function of 
temperature for several bottom quark momenta in the hadron 
gas. 

systems, confirming the prediction of the nonrclativistic 
kinetic theory. These findings constitute an additional 
consistency test of our calculations. 

In the following we study the role of militarization in 
obtaining realistic amplitudes for B meson scattering off 
the light gas and its impact on the transport coefficients. 
We deem this is an important point in order to under- 
stand the wide range of results for transport coefficients 
that can be found in literature within similar approaches 
based on effective theories of heavy meson interactions. 
In Fig. [S] we depict the Bit total cross section in three 
different schemes: using NLO (perturbative) chiral am- 
plitudes, unitarized amplitudes as described in Sec. lIIII B. 
and NLO amplitudes evaluated at threshold energy (scat- 
tering lengths). We observe that the perturbative cross 
section grows monotonically with energy, as expected 
from the chiral expansion. The cross section within the 
threshold approximation is essentially flat over the full 
energy range. The unitarized cross section, in contrast, 
peaks at the resonance region, dominating over the other 
two schemes, whereas later it decreases at higher energies 
as expected from phase space considerations. The corre- 
sponding drag coefficient with p = 0.1 GeV is shown in 



FIG. 8: Upper panel: Bit isospin averaged total cross sec- 
tion for different choices of the scattering amplitudes: uni- 
tarized, perturbative, and scattering lengths. Lower panel: 
Corresponding drag coefficient in static limit as a function of 
temperature. 



the lower panel of Fig. [8] The perturbative scheme leads 
to an unrealistic temperature behavior of F, which is 
clearly tied to the unphysical high-energy behavior of the 
cross section. The scattering length approach does not 
suffer from this artifact at high energies. However, it pro- 
vides a rather smaller diffusion coefficient (factor 1.5 — 2) 
over the whole temperature range as compared to the 
unitarized scheme, since it misses the s-channel enhance- 
ment of the interaction due to the presence of resonances. 
The unitarized scheme leads to the most realistic result 
in our opinion, accounting for the phenomenology of the 
heavy-meson interaction in view of the current knowl- 
edge of the B-meson spectrum, and with a controlled 
high-energy behavior. 

Finally, a comment on the Einstein relation within our 
calculation is in order. In the evaluation of F, Tq and 
Ti we have not made use of the Einstein relation to ob- 
tain the coefficients at low momentum, but performed 
the calculation explicitly as in Eqs. © and ([3]). A perti- 
nent consistency test can be done by obtaining the diffu- 
sion coefficient from the drag coefficient (or viceversa) in 
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static limit from the Einstein relation. Numerically, the 
integrations in Eq. ^ have to be performed with a suit- 
able pion momentum cutoff which ensures convergence. 
One expects the gas distribution function to supress the 
integrand at large pion momentum beyond the chosen 
cutoff so that high-energy contributions are negligible. 
We observe, for instance, that using a cutoff for the pion 
momentum of about 4t:F = 1.2 GcV makes our results 
sensitive at very high temperatures, where the Einstein 
relation is not well fulfilled. On the other hand, one can- 
not increase arbitrarily this cutoff to achieve convergence, 
as this implies that the scattering amplitudes have to be 
evaluated at energies which escape the expected validity 
of the effective theory. Studying how much the static 
coefficients deviate from fulfilling the Einstein relation 
as a function of the momentum cutoff can be considered 
as an estimate of the systematic error in our computa- 
tion of the transport coefficients. An example of this is 
provided in Fig. [§] In the upper panel we plot the func- 
tion F at p = 0.1 GeV (very close to the static limit) 
as obtained from Eqs. (|2I3[) with a momentum cutoff of 
1.2 GcV, together with the determination from To, Ti us- 
ing the Einstein relation. The differences are due to the 
finite pion momentum cutoff. In the lower panel the same 
curves are shown with a cutoff of 3 GeV. This higher cut- 
off ensures convergence of the transport integrals at all 
temperatures and the Einstein relation is well satisfied. 
It is worth mentioning that the scattering length scheme 
discussed above ensures a much quicker convergence of 
the transport integrals, at a cost, however, of implement- 
ing an unrealistic energy dependence of the heavy-light 
meson cross sections and missing the phcnomenological 
information from resonance-enhanced diffusion. 



V. EFFECT OF KAONS AND 77 MESONS 
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FIG. 9: Upper panel: Drag force at p = 0.1 GeV as evaluated 
from Eq. ([2]) and (|30[) . The Einstein relation is not satisfied 
at high temperatures due to the momentum cutoff coming 
from the effectve theory (y^s < 6.5 GeV). Lower panel: Same 
as above but with a larger cutoff momentum. 



Analogously to our previous work, so far we have stud- 
ied the interaction of the bottomed mesons with pions, 
as they are most abundant at heavy-ion collision tem- 
peratures. However, the mesonic gas is also populated 
by kaons and rj mesons, with whom the heavy mesons 
can also interact. In Sec. IIIII we only discussed scatter- 
ing channels involving pions, although we have actually 
extended the computation to the other members of the 
SU(3) octet. With this information it is straightforward 
to implement the contribution from kaons and 77's to the 
transport coefficients. We also added these states for the 
sake of comparison with other references, particularly Q. 

The changes in our calculation are minimal and they 
affect to the collision rate w defined in Eq. ([3]), which 
now reads 

W = Wyr + W K + W-£ + w v . (31) 

We note that the contribution to the transport coeffi- 
cients from different species in the gas is always additive, 
so we expect a moderate increase with the inclusion of 
these new states. 



In Fig. [10] we plot the Tq coefficient of B mesons with 
a momentum p = 0.1 GeV, considering the mesonic gas 
constituted by different particles. We see that the most 
relevant contribution comes from the pion gas, as ex- 
pected. Performing a decomposition of individual hadron 
contributions to To at T = 150 MeV we observe that pi- 
ons provide almost 90% of the total, while the next con- 
tribution is provided by kaons and (mostly) anti-kaons. 
The contribution of r/ mesons is almost negligible, in 
agreement with the fact that the Brj interaction vanishes 
at LO. 

In Fig. [TT] we compare the value of the F and Tq co- 
efficients at p = 0.1 GcV with those of Rcf. jfij]. It can 
be noticed that our results are larger than those of |(|, 
where no unitarization is performed. Tq , though, exhibits 
a smoother growth with temperature in our case, actually 
being overtaken beyond T ~ 150 MeV. We believe that 
the reason of this discrepancy is the use, in Q, of scat- 
tering lengths (valid in principle for very low energies) in 
the whole energy and temperature range. As discussed 
above, this will certainty bring an underestimation of the 
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FIG. 10: To coefficient of B meson with momentum 0.1 GeV, 
considering the mesonic gas constituted by different particles. 



cross section in the region where resonances take over, 
producing lower transport coefficients. In contrast, at 
very high energies, the use of a constant cross section 
may lead to the opposite effect, explaining the behavior 
of r in the high-temperature region from Ref. @ . 

Another point deserving attention is the relaxation of 
bottom quarks in the present situation, for sake of com- 
parison with that obtained in Section HVt which has been 
calculated for a pure pion gas. Taking the meson gas at 
a temperature of 150 MeV, the relaxation length of bot- 
tomed mesons travelling with 1 GeV momentum in the 
full meson gas is 



\ B (T = 150 McV,p = 1 GeV) 



1 

F 



1 



0.012 



fm ~ 83 fm , 



(32) 

We find that the value for Xb in Eq. (|32[) is reduced 
with respect to that estimated in Eq. (|27|) . since in the 
present case the inclusion of other contributions to the 
meson gas yields an increase of the transport coefficients. 
However, this thermal relaxation time continues greater 
than the lifetime of hadron gas, allowing us to consider 
that indeed heavy quarks are carriers of information of 
the phase transition upon exiting the hadron gas. 

For completeness, we also account in Figs. Q21 and Q2] 
for the evolution of the transport coefficients and the ex- 
tracted bottom relaxation length with the heavy-meson 
momentum, in the full meson gas, at a temperature of 
150 MeV. 

In the present approach we can also estimate the en- 
ergy and momentum loss per unit length of a bottomed 
meson travelling in the meson gas from the classical in- 
terpretation of the F coefficient as a drag force. One 
has 



dE/dx = —Fp, and dp/dx=—FE. 



(33) 



in terms of the energy and momentum of the bottomed 
meson. These quantities are depicted in Fig. [H] at T = 
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FIG. 11: F and To coefficient at p = 0.1 GeV in comparison 
with the results of \&. 



150 MeV for the SU(3) meson gas. Thus, a reference 
bottom quark in a B or B* mesonic state with a typical 
momentum of 1 GeV with respect of the rest frame of 
the surrounding medium, loses about 70 MeV per Fermi 
as it propagates in the hadron gas. This effect should 
be taken into account to correctly analyze heavy-meson 
distributions, suggesting that the hadronic phase also has 
to be treated in hydrodynamical simulations regarding 
open heavy flavor probes. 

Finally, we calculate the spatial difusion coefficient, 
D x , which can be related to the diffusion coefficient in 
static limit as D x = T 2 /T (cf. j l] . Appendix B for de- 
tails) . It is shown in Fig. [T5] together with our previ- 
ous result for charmed mesons, below the crossover, and 
with the estimations by Rapp and van Hees for c- and 
fe-quark transport properties above the critical temper- 
ature fl8l - f20l ] . The conclusion observed for the charm 
sector is here confirmed for bottom: the minimum relax- 
ation time for heavy flavor seems to take place around the 
crossover, where one expects strongest (and long-range) 
interactions. Our present work reinforces the use of heavy 
flavor as probes of the QCD phase transition by reduc- 
ing theoretical uncertainties regarding the dynamics in 
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the hadronic phase. 



VI. SUMMARY AND CONCLUSIONS 

In this work the drag and diffusion coefficients of bot- 
tomed mesons in a thermal gas of pions, kaons and rj 
mesons have been evaluated in a Fokkcr-Plank transport 
approach. The dynamics of the interaction of B and 
B* with the light mesons has been modelled employing 
a unitarized version of heavy-meson chiral perturbation 
theory within the constraints of heavy quark symmetry. 
The relevant scattering amplitudes have been calculated 
at next-to-leading order in the chiral expansion and lead- 
ing order in the heavy quark limit, and the free parame- 
ters have been constrained by available data on the B^ 
spectrum and heavy-quark considerations. 

We have observed a sizable temperature dependence 
of the transport coefficients, indicating that both drag 
and diffusion arc more efficient in the hotter stages of 
the hadronic phase in a heavy-ion collision. 

Regarding momentum dependence, it turns out that 
F and Tq only change mildly whereas transverse coeffi- 
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FIG. 13: Bottom relaxation length in the hadronic phase at 
p = 0.1 GeV as a function of the temperature. 



cient, Ti, is particularly sensitive to the momentum of 
the heavy meson. The latter may lead to observable con- 
sequences in the analysis of anisotropic observables such 
as the elliptic flow. 

In addition, a detailed comparison has been accom- 
plished between our outcomes and the other ones in liter- 
ature obtained via different approaches. Our statements 
rely on a good control of the elementary B-meson inter- 
action with the light meson octet over a wide range of en- 
ergies, accounting for dynamical generation of resonances 
in attractive channels. We have shown that preserving 
unitarity in the scattering amplitudes plays an esential 
role to provide realistic estimates of the transport coeffi- 
cients at high temperatures. 

We have performed several consistency tests of our re- 
sults such as studying the scaling of transport coefficients 
with the heavy-quark mass and verifying the Einstein re- 
lation in the static limit, which were satisfactory. 

Moreover, the individual contribution of pions, kaons 
and r\ mesons to the transport coefficients have been es- 
timated. The most relevant contribution is lead by the 
pion gas, as expected, with the next-to-leading contribu- 
tion being provided by kaons and anti-kaons. 

Also, some estimations of relevant quantities have been 
done directly from the transport coefficients. One exam- 
ple is the bottom relaxation length, which is about 83 
fm for bottomed mesons with momenta of 1 GeV and for 
meson gas at temperature of 150 MeV. It allows us to in- 
fer that bottomed mesons barely relax during the lifetime 
of the hadron gas, unlike charm mesons that, while not 
relaxing completely, they loose a great deal of memory of 
the initial state. In this sense, the bottomed mesons con- 
stitute an optimal system to characterize the early stages 
of a rclativistic heavy-ion collision. Observables like the 
nuclear suppresion factor or the elliptic flow can provide 
clear indications of the system properties and evolution 
after the nuclear collision. 

Furthermore, the computation of the spatial diffusion 
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tomed mesons at a temperature of 150 MeV in the full meson 
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coefficient for the present case is in agreement with the 
idea that the relaxation time for heavy quarks has a min- 
imum around the crossover to the quark-gluon plasma, 
the expected place of strongest interactions. 

Another interesting quantity evaluated has been the 
momentum loss per unit length, which is 70 MeV per 
Fermi for a bottomed meson with momentum of 1 GeV 
propagating in the hadron gas. This result indicates that 
to correctly analyze the heavy-meson distributions, the 
effect of loss of energy and momentum must be taken into 
account. 

Hence, the findings of the present work above dis- 
cussed reinforce the role of heavy mesons as probes of the 
strongly interacting matter phase transition, paving the 
way to a better understanding of heavy-flavor dynamics 
and transport properties below the crossover. 
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Appendix A: Subtraction constant 

In paralcll with the work [3| one would like to pro- 
vide a reasonable description of meson-meson scattering 
in such a way that the pertinent resonances in the dif- 
ferent channels are generated. In particular, we would 
like to reproduce the Bq and B\ resonances in the Bit 
and B*tt channels, respectively. However, the masses 
and widths of these resonances are not yet experimen- 
tally known. 

Starting at LO, there are no free low energy constants 
and the perturbative amplitude V in the B*tt channel is 
fixed. In the unitarized amplitude -the one which gener- 
ates a pole in the amplitude- one introduces a subtrac- 
tion constant a(jx). This constant will be fixed by match- 
ing the function G(s) in dimensional regularization to the 
one with cutoff regularization scheme. 

The loop function G(s) in dimensional regularization 
is given in Eq. (fTi?|) . The regularization scale /j, is chosen 
at a natural scale of 1 GeV [25| In a cutoff regularization 
one has [2lJ 
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FIG. 16: Re G(s) and lm G(s) from Eq. (H) (dimensional 
regularization) and Eq. (|A1() (cutoff regularization) . We 
match the loop function in the two regularization schemes 
at •jsth = niB + rn^. A higher energies G A (s) becomes non- 
physical and diverges around y/Soo = 6 GeV. 



Note that this expression is valid below and above 
threshold. The equivalent cutoff in momentum, A, is ob- 
tained by matching both expressions where we use the 
value of a(fi) in the Dir sector It is well known 
that Ga presents a spurious divergence at some above 
threshold, where its determination is unreliable. For this 
reason, one typically compares the two loop functions at 
threshold s t h = {tub + m^) 2 (see Fig. H"6|). 



Combining all these requirements we obtain 
a(l GeV) = —3.47, corresponding to a reasonable 
cutoff momentum of about 1 GeV. Within this scheme 
we obtain a resonance around 5580 MeV at LO in the 
7 = 1/2 channel, about 100 MeV below the reference 
value of the mass of the B\ (heavy quark spin 3/2). 
At NLO the two free low energy constants h 3 and /15 
can modify the pole position and width. Varying these 
parameters within the constraints of HQ Symmetry, 
a maximal value of 5587 MeV can be found for the 
pole position, with a width around 245 MeV. Although 
higher values of the resonance mass can be forced by 
tuning the subtraction constant, in order to keep the 
equivalent cutoff of natural scale we content ourselves 
with a .Bi pole mass of 5587 MeV and 5534 MeV for the 
Bq resonance (with a width of 210 MeV). Our results 
agree with previous studies using similar unitarization 
methods. They are summarized in Table [TTTl 
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